Fine structure of the local pseudogap and Fano effect 
for superconducting electrons near a zigzag graphene edge 
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Motivated by recent scanning tunneling experiments on zigzag-terminated graphene this paper 
investigates an interplay of evanescent and extended quasiparticle states in the local density of states 
(LDOS) near a zigzag edge using the Green's function of the Dirac equation. A model system is 
considered where the local electronic structure near the edge influences transport of both normal 
and superconducting electrons via a Fano resonance. In particular, the temperature enhancement 
of the critical Josephson current and — tt transitions are predicted. 

PACS numbers: 73.23.Ad,74.50.+r, 74.78. Na 

Introduction.- Experimental evidence^'^ for massless 
Dirac-like quasiparticles in graphene - a carbon mono- 
layer with the hexagonal structure - has stimulated vig- 
orous interest in electronic properties of this system (e.g. 
Refs. iiUlQSiliOini)- The unit cell of graphene con- 
tains two atoms each belonging to a triangular sub- 
lattice, and the low-energy states are described by a 
two-dimensional Dirac equation where the role of spin 
is assumed by the sublattice degree of freedom (pseu- 
dospin)i^ii^. Similar to relativistic spin-half particles in 
two dimensions, the graphene bulk density of states has 
a linear pseudogap^ around zero energy E = 0. Natural 
boundaries can however give rise to additional spectral 
branches such as the low-energy edge states^ii^. They 
are localized near a zigzag-shaped edge, whose outermost 
sites all belong to the same sublattice [Fig. UK a)], and 
originate from the effective pseudospin "polarization" 
due to vanishing of one of the pseudospinor components 
as required by particle conservation^. Recent scanning 
tunneling experiment o^^'^^ report a singular enhance- 
ment of the LDOS near zigzag boundaries attributed to 
the edge states. 

The measurementsiSiii also revealed another peculiar- 
ity of the energy dependence of the LDOS - a fine os- 
cillatory structure superimposed on the pseudogap with 
the amplitude enhanced at larger energies^''. The origin 
of this behavior is still unaccounted for, although sub- 
sequent publications have studied the graphene LDOS, 
e.g. numerical simulations of Ref. [l^ found damped spa- 
tial oscillations of the LDOS. The present study indends 
to show that both findings are consistent with the picture 
of interfering Dirac electron waves near a zigzag edge. To 
demonstrate this point, the one-particle Green's function 
of the Dirac equation was calculated for clean graphene 
with a zigzag edge described by the boundary condition 
of Ref. [1. Then, the following expression for the LDOS 
v{E,d), as a function of energy E and distance d from 
the edge, was obtained 



components given by the Bessel functions jQ{2Ed/hv) 
and Ji(2\E\d/hv) are due to interfering waves formed 
of the states belonging to the Dirac spectrum {v and 
h are the electron velocity and Planck's constant). For 
2\E\d/hv ^ 1 the amplitude of the oscillations is pro- 
portional to y^\E\/d [see also Figs. [Ijb) and (c)], which 
qualitatively agrees with both the experiment— and nu- 
merical simulationsi^. 

Another issue this study focuses on is the connec- 
tion between the local electronic structure of zigzag- 
terminated graphene and Fano scatteringi^i. Unlike ear- 
lier works [e.g. Ref.H^I where the Fano effect was due to 
resonant flux states in finite-size ribbons, here the Fano 
resonance is studied in a nanowire side-coupled to half- 
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Here the delta-functional term results from the disper- 
sionless zero-energy edge state whereas the oscillating 



FIG. 1: (Color online) (a) Schematic view of a zigzag 
graphene ribbon terminated by atomic lines belonging to dif- 
ferent sublattices, conventionally denoted as A and B. (b) 
LDOS vs. energy at different distances from the edge: (A) 
d = 4a, (B) d = 10a, (C) d = 30a, (D) d = 10000a; where 
a = 0.246 nm is graphene's lattice constant, Va ~ 1/Snhva 
and Ea ~ hv/2a. (c) LDOS vs. distance from the edge for 
different energies: (A) E = 0.15^, (B) E = O.ZEa, (C) 
E — 0.5Ea- The delta function in Eq. is approximated by 
a Lorenzian o/tv{E^ + o^) with o = 0.03i?a. 
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infinite graphene, by analogy with similar quantum-dot 
structures^^i^i^i, and originates from a single disper- 
sionless edge state. Also, unlike Ref. [l^, the main focus 
here is on the transport of correlated electrons in Joseph- 
son nanowires. The Fano effect is predicted to cause quite 
unusual behaviors of the critical current such as enhance- 
ment by temperature and, under certain conditions, — tt 
transitions similar to those in ferromagnetic Josephson 



junctions 
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In the context of the Josephson 



effect in graphene nanostructures (e.g. Refs. lslHs^ ) these 
issues have not yet been addressed. 

Green's function of a zigzag ribbon.- Assuming 
no scattering between the two valleys, K and K\ of 
graphene's Brillouin zone^^, one only needs to calcu- 
late the Green's function in one of them, e.g. K, where 
the Dirac equation reads [uqE -t- ihv{axdx + aydy)\G = 
crQ5{x — x')S{y — y'). Here the (retarded) Green's func- 
tion matrix Gjk with j, k — A, B, Pauli ax,y and unity 
CTq matrices all act in pseudospin space. It suffices to 
solve the pair of equations for Gaa and Gba- After 
expanding in plane waves e^^^ , the equations for the 
Fourier components are Gg^lk = {fiv/E){k + dy)GAA\k 
and [dl - q^]GAA\k = {Elh^v^)5{y - y') with = 
— (E/Hv)'^. The solution can be sought in the form 
GAA\k{y, y') = a{y')e-^y + b{y')eiy - Ee-'^^y-y'^ /2h^v^q, 
where the last term is the Green's function of an un- 
bounded system, and the coefficients a(i/') and b{y') are 
to be found from the boundary conditions^ GBA\k\y=o = 
{k + dy)GAA\k\y=o = and GAA\k\v=w = 0. This yields 
the following result 

GAA\kiy,y')=E/2h\^qx 

k[coshq{w — |y — y'\) — coshq{w — y ~ y')] 



q[smhq{i 



q cosh qw — k sinh qw 
~\y~ y'\) + sinhg(w -y- y')] 
q cosh qw — k sinh qw 



(2) 



The poles of GAA\k, given by the equation q = ktanh qw 
(cf. Ref. d), determine the excitation spectrum. As 
knowni^ii^, it has an almost flat branch merging with 
the Fermi level E — corresponding to a state exponen- 
tially decaying from the edge into the interior. This can 
be easily seen from Eq. ([2]) in the limit w — > oo: 



GAA\k{y,y') 



Ec-i\y-y'\ (g + A:)2e-9fe+y') 



(3) 



The pole E — describes a dispersionless edge state 
existing for /c > 0. From Eq. ([3]) an exact position 
representation for the Green's function GAA^xy, xy') = 
J'^^dkGAA\k{yjy')/{'^'^) can be obtained as 



GAA{xy,xy') = 



EYoikE\y-y'\)-i\E\MkE\y~y'\) 



{2hvY 

EYojkEiy + y')) - i\E\JoikEiy + y')) 
{2hvY 

2EY^{kE{y + y')) - 2i\E\J^{kE{y + y')) 
{2hvYkE{y + y') 



(4) 




FIG. 2: (a) Schematic view of a zigzag-edge grapliene ribbon 
witli a side tunnel contact to a nanowire connecting electron 
reservoirs 1 and 2. The contact is assumed point-like, i.e. 
its size is much bigger than the interatomic distances, but 
smaller than the electronic mean free paths in both wire and 
graphene. (b) Zero-temperature conductance g (in units of 
/-kK) vs. bias voltage V for spin-degenerate ih — 0) and 
spin-split (h — ST) edge state in graphene. 



where Jniz) and l^(z) (n = 0,1) are, respec- 
tively, the Bessel and Neumann functions, and ks = 
Vs^/hv. The LDOS [Eq. (P)] is obtained via i^{E,d) = 
— {l/Tr)lmGAA{xy,xy')\y=y'=d, taking into account the 
pole of the function Yi. We note that the interference 
of the waves incident at and reflected from the edge with 
small momenta |fc| < \E\/hv produces spatial oscillations 
of the LDOS with the period much bigger than the lat- 
tice constant a [Fig. [TI^c)], i.e. well within the scanning 
tunneling microscop resolution. It is also instructive to 
examine Eq. (j4|) near the edge where it assumes a uni- 
versal form. 



GAA{xy,xy') w l/TrE{y + y') 



l\2 
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0, 



(5) 



independent of material parameters. To regulate the 
divergence a.t y — y' = 0, due to the effective contin- 
uum description, it is convenient to introduce the cutoff 
GAAixO, xO) w l/AnEdl with 4 ^ a- 

Fano scattering off a zigzag edge.- The behavior 
of the Green's function near the edge can have a direct 
impact on charge transport. Let us consider a quasi-one- 
dimcnsional wire (with conventional quasiparticlc spec- 
trum) coupled in parallel to a zigzag graphene edge via a 
point-like tunnel barrier [Fig. [5]Ja)]. The contact is mod- 
elled by a real-space tunneling Hamiltonian of the form 
Ht = Vi(0)(rV'A(ro) +rV;i(ro)) +h.c., ro = (0,0), 
where the electron operator in the wire -ipiiO) at the 
contact point a; = is coupled to those in graphene on 
sublattice A in both valleys K, tpAiJ^o) and K', ?/'^(ro) 
with the matrix elements T and T'. To describe electron 
scattering caused by the contact, I use the equations-of- 
motion method and calculate the retarded Green's func- 
tion in the wire 



i-gL°)(o,o)e 

Here Gw\x,x') — e*''"™!^^^ ^/ikvp is the Green's func- 
tion in the absence of tunneling (kw ~ fci? -I- E/hvp with 
vp and kp being the Fermi velocity and wave number 
in the channel), S = \T\^GAA{ro,ro) + \T'\^G'^^{ro,ro) 



FIG. 3: Critical current vs. temperature for a spin- 
degenerate edge state: h — 0, j = O.IEl, A = WEl- The 
current is normalized to the value IciTmin) where Tmin ~ 
0.05-EL/fes is the lowest temperature for which the condition 
7 < irkBT ^ A of weak proximity effect still holds. 

is the tunneling self-energy, and G^^(ro,ro) is the 
Green's function in valley K' coinciding with Gaa{t^q, to) 
[Eq. ([5])]. The transmission amplitude''^ between the 
reservoirs is t = ihypGw (f = e''^^^E/{E + iT), 

where T = (|Tp -I- \T'\'^) / AndlKvp determines the tun- 
neling rate T/fi between the systems. The Fano-like 
transmission antiresonance at -E = manifests complete 
backscattering of an electron wave incoming from one 
of the reservoirs. It is due to destructive interference be- 
tween the electron wave directly transmitted through the 
wire (without tunneling) and the wave transmitted via 
tunneling through the graphene edge state whose energy 
is pinned to the Fermi level in the wire. It is straightfor- 
ward to generalize the analysis to a spin-split edge state 
with energies =f/i for spin projections a — ±1/2. In this 
case, we have 

t„{E) = e'''^^{E + 2ah)/{E + 2ah + iV). (7) 

Figure [DJb) shows the voltage dependence of 
the zero-temperature Landauer conductance 
g{V) = eV(27r?i)Ea=±i/2l^a(e^)l'- For ;^ ^ 
the conductance dip is split due to the spin-filtering 
effect discussed earlie r'^'^i^^ in the context of possible 
applications in spintronics^. 

Fano effect in a Josephson junction.- Let us fi- 
nally discuss the case of superconducting reservoirs sup- 
porting an equilibrium Josephson current. The Joseph- 
son coupling is maintained due to the Andreev process^ 
whereby an electron is retro-reflected as a Fermi-sea hole 
from one of the superconductors with the subsequent 
hole-to-electron conversion in the other one. Such an 
Andreev reflection circle facilitates a Cooper pair trans- 
fer between the superconductors. Since both electron 
and hole also experience normal scattering inside the 
junction, the transmission antiresonance is expected to 
strongly influence the Josephson current. It is conve- 
nient to use the approach of Refs. .38.i39., relating the su- 
percurrent to the scattering amplitudes via a sum over 
the Matsubara frequencies w„ = {2n-\- l)TrkBT as follows 

I, = -AekBT/h J2 B.liE)t„{E)t*_^{-E)\E=^u^, (8) 

n>0,a 



FIG. 4: Critical current vs. spin-splitting energy: T = 
O.lSi/fcs, 7 = O.IEl, A = WEl. 

where t*ia{—E) is the hole transmission amplitude cor- 
responding to the time-reversed counterpart of the elec- 
tron Hamiltonian^Sii^, and aa{E) is the Andreev reflec- 
tion amplitude at the point contacts to superconductors 
1 and 2. Equation ([8]) is applicable for arbitrary ta{E) 
as long as a^ is small enough so that one can neglect 
higher order Andreev processes. In this case the Joseph- 
son current-phase relation is sinusoidal with Ic [Eq. ([5])] 
being the critical value of the current. 

In contacts to conventional Bardeen-Cooper-Schrieffer 
(BSC) superconductors, the Andreev process is well 
described by the scattering model of Ref. |4l|. How- 
ever, in many practical cases superconducting contacts 
to low-dimensional systems can hardly be regarded as 
BCS-like ones. Proximity-effect contacts to semiconduc- 
tor nanowires^2ii^ and carbon nanotubes^i^ are impor- 
tant examples of such a situation. In this case a thin 
normal-metal layer is inserted between the superconduc- 
tor and the wire to ensure a good electrical contact. 
In proximity-effect point contacts the Andreev scatter- 
ing amplitude can be expressed in terms of the qua- 
siclassical condensate J-a{oJn) and quasiparticle Gai^n) 
Green's functions of the normal layer a,s^^^ aa{ujn) = 
iJ-a/{l + Qa). I will adopt this approach and make 
use of McMillan's expressions^^i^ for the Green's func- 
tions: Ta = A„/y^^TAf, = {uJn/An)J^a, and 
A„ — 7 A/ (7 + ^cj2 + A^), where A is the supercon- 
ductor's pairing energy and 7 is McMillan's parameter 
controlling the strength of the proximity effect in the nor- 
mal layer and, hence, the Andreev reflection amplitude 
&ai^n) = i^n/i^n + •\/a;^ + A^ ). For a weak proximity 
effect with 7 < TrfcsT <C A, the amplitude a^ is small^ 
and equation ([5]) assumes the form 

' n K + V^J^Ta^]^ '[/. + z(.„ + r)p' 

where the exponential factor results from the dynami- 
cal phase 2EL/hvF accumulated in the Andreev circle, 
introducing the Thouless energy El = hvF/2L. 

Figure O shows the temperature dependence of Ic for 
the spin-degenerate case (/i = 0). In the absence of tun- 
neling (r = 0) it is just a monotonic exponential de- 
crease. However, for F 7^ the interplay of the trans- 
mission antiresonance and the exponential suppression 
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FIG. 5: Critical current vs. temperature for a spin-polarized 
edge state: h = 0.5El, 7 = O.IEl, A = IOEl- 

gives rise to a maximum at finite T. Spin splitting of 
the graphene edge state lifts the Fano resonance condi- 
tion E — ior both electron and Andreev reflected hole. 
Therefore, for relatively weak tunneling coupling, when 
■KksT < r < the critical current increases with h 
[Fig.S], a behavior quite unusual for Josephson junctions. 
Surprisingly, for stronger tunneling coupling (F > E^) 
the function Ic{h) becomes nonmonotonic with a rather 



broad region hi < < /i2 ~ T where Ic is negative. The 
lower boundary hi « TThgT is set by the temperature and 
is much smaller than all El, F and A. The supercurrent 
reversal is a consequence of the spin-dependent phases 
acquired by both electron and hole due to scattering off 
the spin-polarized graphene edge, with negative values 
of Ic implying a built-in 7r-phase difference in the ground 
state of a Josephson junction^^ as opposed to 0-phase dif- 
ference for Ic > 0. The — tt transition can be driven by 
temperature as shown in Fig. [5l Such a tt state is known 
to occur in ferromagnetic junctions where the conden- 
sate function oscillates in space2^iS£i^i2^ (see, also recent 
reviews [29ll30l ) . The author is not aware of any earlier 
work predicting — tt transitions due to spin-dependent 
Fano scattering. The main condition for this mechanism 
to work, i.e. \h\ > nkBT, can be met in the millikelvin 
region at modest external magnetic fields. 
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